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By using 14 different versions and parametrizations of a proximity potential and two new versions 
of the potential proposed in this paper, we perform a comparative study of fusion barriers by 

studying 26 symmetric reactions. The mass asymmetry tja = ^ A2+A* ) > however, is very large. 
Our detailed investigation reveals that most of the proximity potentials reproduce experimental 
data within ±8% on the average. A comparison of fusion cross sections indicates that Bass 80, 
AW 95, and Denisov DP potentials have a better edge than other potentials. We also propose new 
versions of the proximity potential as well as Denisov parametrized potential. These new versions 
improve the agreement with the data. 
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INTRODUCTION 



The study of fusion barriers and fusion cross sections has received renewed attention in recent years [lH5(. This 
renewed interest is caused by the efforts of low and intermediate energies in order to understand the fusion mechanism 
■ and, subsequently, the nucleus-nucleus interactions in nuclear physics. This is further boosted by the availability of 
radioactive- ion beams involving various reactions Q . 

In recent years, large numbers of models depending on the vast variety of assumptions have been proposed [HEHHS - 
^1 ' H3] • One set of such theoretical models is based on the microscopicpicture in which one starts from two- and three- 
body effective interactions and calculates the ion- ion potential (3, |j, H| - Another class of the models takes the gross 
macroscopic picture into account [III EH- Because of the recent precise measurements of the fusion cross sections, 
the job of theoretical models has become much more complex. Among different theoretical models, the proximity 
potential enjoys special status [1, 0|- All proximity potentials are based on the proximity force theorem, according 
to which, the nuclear part of the interaction potential can be written as the product of a factor depending on the 
mean curvature of the interaction surface and a universal function (depending on the separation distance) and is 
independent of the masses of colliding nuclei. This concept did introduce a great amount of simplification in nuclear 
potential studies [U, H| ■ Several refinements and modifications have been proposed in the recent past concerning over 
the original proximity potential to remove the gray part of the potential 6] . This observation demands a careful and 
systematic study of the heavy-ion fusion process using all such potentials that are remodeled and parametrized within 
the proximity concept. 

Here, we will concentrate on the symmetric colliding pairs with symmetry parameter A s =( T ) = 0.0 only. Here, 



■A 

\S 1 N and Z belong to the combined neutron and proton content of the reaction. This comparison, which covers a wide 
spectrum, will give us a unique possibility to compare various modeled proximity potentials. We also plan to modify 
the proximity potential and the potential by Denisov and we will show that these new versions improve the agreement 
with the experimental data. Section |H] describes the Formalism in brief, Section Mil depicts the results and Summary 
is presented in Section ITVl 

II. FORMALISM: 

According to the original version of the proximity potential 1977 [7|j, the interaction potential VXr(r) between two 
surfaces can be written as 

yProx 77 ( r ) = 47r7 57j$ (r - Ci - C 2 /6) MeV, (1) 
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where 7 is the surface energy coefficient. The mean curvature radius, R in Eq. (Q]) is written as: 

R = CiC 2 /(d + C 2 ); C t =Ri[l- (b/R t ) 2 + •■•], (2) 
and Rj, the effective sharp radius, reads as; 

Ri = 1.28A 1 / 3 - 0.76 + 0.8A~ 1/3 fm (i = 1, 2). (3) 

This model is referred to as Prox 77, and the corresponding potential is referred to as V^ rox 77 (r). 

Later, Reisdorf [8| modified the preceding potential with a different 7 value. This is labeled as Prox 88. 

The latest version of the proximity potential by Myers and Swiatecki (labeled here as Prox 00) is given in Ref. [f| . 
This Prox 00 uses an experimental value and/or an old formula for the radius. 

Recently, Royer and Rousseau [2] gave a more precise radius formula 

Ri = aA] /3 {1 + 0/Ai - SA sl } fm (i = 1,2), (4) 

where a, /?, and S are constants that have values of 1.2332, 2.348443 and 0.151541, respectively. This formula is 
obtained by analyzing as many as 2027 masses with N, Z > 8 and a mass uncertainty < 150 keV. We implement this 
radius formula in Prox 00 and label it as Prox 00DP. 

Based on the proximity concept, many other potentials have also been shown in the literature. We will use the 
potentials by Bass 1973 (labeled here as Bass 73) [nj], similarly, Bass 1977 (Bass 77) [12], Bass 1980 (Bass 80) [j, 
Christensen and Winther 1976 (CW 76) [3, Broglia and Winther 1991 (BW 91) 0, Aage Winther (AW 95) [H|, 
Ngo 1975 (Ngo 75) 0], Ngo 1980 (Ngo 80) [ij], and Denisov [|. 

The potential by Denisov @ is also modified here to include the previous more precise radius formula [Eq. @] in 
its parametrization. This is labeled as Denisov DP. 

The exact potential based on the Skyrme energy density formalism [3j (labeled as EDF Exact) along with its 
parametrized form (labeled as EDF Par) [r| will also be used for comparison. Note that Skyrme forces are also 
widely used in intermediate energies p~7j | . 



III. RESULTS AND DISCUSSIONS 



In Fig. 1, we display the nuclear part of the interaction potential Vjv(MeV) as a function of internuclear distance 
r (in femtometers) using some of the previously listed versions of the proximity potentials for the reactions of 16 0+ 16 
and 40 Ca+ 40 Ca. In Fig. 1(a), we display three versions of the potentials by Bass, whereas in Fig. 1(b) we deals 
with Ngo parametrizations. The potentials of Winther and collaborators are displayed in Fig. 1(c), followed by four 
versions of a proximity potential displayed in Fig. 1(d). 

From the figure, we see that different versions of the Bass potentials as well as the CW 76 do not have a repulsive 
core at shorter distance. On the other hand, the BW 91 and AW 95 potentials, follow the Woods-Saxon-type 
distributions. All other potentials have acceptable shape: that is, attractive at long distances followed by repulsive 
at shorter distances. Interestingly, the deepest potential caused by the proximity potential Prox 88 (=72 MeV). In 
other words, one see a huge difference in the potentials obtained from different versions of proximity-based formalisms 
even within the same model. For instance, different versions of the proximity potential differ by as much as 17 MeV 
for the reaction of 40 Ca+ 40 Ca. Since the fusion process is a low-density phenomenon happening at the outer surface, 
the difference in the inner part of the potential may not be so important. By adding the Coulomb potential to the 
nuclear part, one can compute the total potential Vr(r) as 

V T (r)=V N (r) + V c (r), (5) 

= Mr) + ^. (6) 
r 

Since the fusion happens at a distance greater than the touching configuration of the colliding pair, the previous form 
of the Coulomb potential is justified. One can extract the barrier height Vg heor and the barrier position Rg leor using 
the following conditions 

dV T (r), n , d 2 V T {r). 

—^±\ r=Rror = 0; and _^Ai| r=fl&w < 0. (7) 

In Fig. 2, we display the total interaction (nuclear + Coulomb) potential Vr(MeV) as a function of internuclear 
distance "r" (in femtometers) . We display the results using different versions of the proximity potentials (Prox 77, 
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Prox 88, Prox 00, and Prox 00DP) and potentials of Winther and collaborators i.e., (CW 76, BW 91, and AW 95). 
From Fig. 2, we note that, although various potentials differ significantly in the interior part, very little dependence is 
visible on the surface part and also on the barrier region. In terms of fusion barrier heights and positions, all versions 
yield nearly the same barrier heights and positions. Of course, the shape and the curvature of the potential differs, 
indicating a different picture for subbarrier fusion cross sections that depends sensitively on the shape of the potential 
one is using. For the reaction of 40 Ca + 40 Ca, no difference is seen between Prox 00 and Prox 00DP. This happens 
because the experimental value of the charge distribution is available for the 40 Ca nucleus, and both models use the 
same experimental values. The results may differ for those nuclei in which experimental values are not available. 

Knowledge of the shape of the potential, as well as the barrier position and the height, allows one to calculate the 
fusion cross section at a microscopic level. To study the fusion cross sections, we will use the model given by Wong 
[l8j . In this formalism, the cross section for complete fusion is given by 



where k — y / 2/iE/H 2 and here [i is the reduced mass. The center-of-mass energy is denoted by E cm . In this formula, 
Imax corresponds to the largest partial wave for which a pocket still exists in the interaction potential, and Tj (E cm ) 
is the energy-dependent barrier penetration factor and is given by, 



Ti (E cm ) = 1 + cxp 



2tt 



\rtheor rp \ 



(9) 



where hw\ is the curvature of the inverted parabola. If we assume that the barrier position and the width are 
independent of /, the fusion cross section reduces to 



a fus (mb) = 



lQ R theor fa 
- X 

^E cm 

In < 1 + exp 



2vr 

hujQ 



(E cm - Vjf eor ) 



(10) 



For E crn >>V^ eor , the preceding formula reduces to a well-known sharp cutoff formula 

rthe 



t in 



cT fus (mb) = 10ttB% 



1 



Earn. 



(11) 



whereas for E cm <<V^ eor , the foregoing formula reduces to 



Of us (mb) 



IQRtheor' 



2E r 



V, 



theor\ 



(12) 



We used Eq. (jT0|) to calculate the fusion cross sections. 

In Fig. 3, we display the theoretical barrier heights V^ heor (MeV) verses the experimental barrier heights Vg Xpt (MeV) 
using all 16 different potentials. The experimental barrier heights V^ xpt are taken directly from the literature di- 
rectly |19l - l25| . The limited numbers of reactions in certain cases are caused by the restrictions posed on different 



potentials [111, (TM^ • The lines are the fits over the points. These fitted equations distinguish tell the deviation from 
the experimental data. Very interestingly, we see that, all models can reproduce the experimental barrier heights for 
symmetric colliding nuclei within ±8% on average. We also notice that, on average, barriers formed using EDF Exact, 
EDF Par, Bass 77, Bass 80, Denisov DP, and the different versions by Winther and collaborators are close to the 
experimental data. On the other hand, barriers formed within Bass 73, Ngo 80 and Prox 77 potentials deviate by ± 
7% from the experimental values. The revised versions by Bass improve the barrier heights drastically. Now Bass 77 
and Bass 80 reproduce the experimental data within 1.5%. A newer version of the Bass potential (Bass 80) shows 
slight improvement over Bass 77. Strangely, the Ngo 80 version deviates more than 7%, whereas its first version, was 
able to reproduce the barrier heights within 2%. For Ngo 75, only two systems fall within its parametrization limits. 

The four versions of the proximity potentials yield an interesting comparison. As pointed out by various authors 0] , 
Prox 77 deviates from the experimental fusion barrier heights by 4% (here it is 6.73%) whereas, the new versions of 
the proximity results are lowered to 0.01% (here, it is 5.34%). Note that in Ref. Q, very old data for the fusion 
barriers were used. In this paper, we have used the latest data. It is worth mentioning that a slight variation in the 
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radius formula (Prox 00 and Prox 00DP), can improve the comparison by nearly 1%. Very interestinglya change in 
the value of the surface energy coefficient 7 (Prox 77 and Prox 88) improves the agreement drastically |26|. Further, 
the advantage of a new proximity potential (Prox 00) can clearly be obtained by just using 7 in Prox 77 as suggested 
by Reisdorf [8J . We do not see a direct advantage of the original Denisov form of the potential over the Skyrme energy 
density model by Puri and Gupta |16| in which a perfect comparison with the experimental data is clearly visible. Its 
new form, Denisov DP, however, shows perfect agreement with the experimental data. From these figures and this 
analysis, it is very clear that different models do not yield very different results. Instead, technical parameters such 
as the surface energy coefficient 7 and the radius, can have a significant impact on the outcome. The implementation 
of the latest radius formula clearly yields better agreement. 

In Fig. 4, we display the fusion barrier positions as a function of experimentally extracted values. We see that no 
trend emerges in this case. This happens because of a great amount of uncertainty in the measurement of the fusion 
barrier positions reported by various authors in various experiments pol . l2lj . 

In Fig. 5, we display the percentage difference of the fusion barrier heights over its experimental values defined as; 

ytheor _ yexpt 

^V B (%) = x 100. (13) 

V B 

We see that on an individual bases, all proximity potentials can reproduce the data within ±10%. The least amount 
of deviation is attained from Bass 77, Bass 80, different Winther potentials (CW 76, BW 91, and AW 95), EDF 
Exact, EDF Par and Denisov DP ( all within ±5%). However, the proximity potential (Prox 77) deviates much more 
compared to its other versions. 

In Fig. 6, we display the fusion cross section Uf us (in millibarns) as a function of center-of-mass energy E cm for 
the reactions of 24 Mg + 28 Si [Fig. 6((a)] and 40 Ca + 40 Ca [Fig. 6(b)]. Here, the latest versions of the proximity 
parametrizations along with the original proximity potential and its modifications are shown for clarity. The exper- 
imental data are taken from the Refs. (2lT - l25l |. As we see, Bass 80 Denisov DP, and AW 95 do a better job for the 
reaction of 24 Mg + 28 Si whereas Prox 77 and Ngo 80 fail to come close to the experimental data. For the reaction of 
40 Ca + 40 Ca, no clear picture emerges. In both cases, the potential of Denisov DP and AW 95 are able to reproduce 
the cross section. 



IV. SUMMARY 

By using as many as 16 versions of the proximity potential derived either from the proximity potential or from 
the parametrized versions in terms of the proximity concept, we carried out a comparative study of fusion barriers 
for symmetric colliding nuclei. For the present study, four versions of the proximity potential, three versions of the 
proximity potential by Bass, three versions of the proximity potentials by Winther and collaborators, two versions of 
the proximity potentials by Ngo and two versions of the proximity potentials by Denisov and EDF each were taken. 
We also proposed new versions of the proximity potential and a proximity potential by Denisov. A detailed study 
reveals that all potentials can reproduce experimental data, on average, within ±8%. However, the comparison of 
fusion cross sections reveals that the Bass 80 Denisov DP, and AW 95 potentials reproduce data better than the other 
potentials. 
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r (fm) 



FIG. 1: (Color online) The nuclear part of the interaction potential, VW(MeV) as a function of internuclear distance r (in 
femtometers) for the reactions of 16 O + 16 O and 40 Ca + 40 Ca using different proximity potentials. 
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FIG. 2: (Color online) The total interaction potential Vr(MeV) as a function of internuclear distance r (in femtometers). Here, 
we display the results obtained with different versions of proximity and Winther potentials only. 
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FIG. 3: Comparison of theoretical fusion barrier heights Vg heor (MeV) using different proximity potentials with experimental 
values Vg Xpt (MeV) 19 25]. The solid lines represent the straight line least squares fit created over different points. 
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FIG. 4: A comparison of theoretical _RB eor (fm) and experimental fusion barrier positions 7?^ p *(fm) [19 - [2a ] using various 
versions of the proximity potential. Solid line represent the straight line least square fit. 
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FIG. 5: The percentage deviation AVb (%) as a function of the product of charges Z1Z2 using different versions of the proximity 
potential. 
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FIG. 6: (Color online)The fusion cross-sections for the reactions of 24 Mg + 28 Si [Fig. 6(a)] and 40 Ca + 40 Ca [Fig. 6(b)] as 
a function of center-of-mass energy E cm . The experimental data are taken from Morsad 1990 [22], Gary 1982 [23j], Tomasi 
1982 [53], Aljuwair 1984 [53], and Barreto 1983 [2^|. For the clarity, only the latest versions of the different proximity potentials 
are shown. 



